FIELD OF THE INVENTION 

The present invention relates to a method referred to as lumping or pseudoisation 
method for describing a fluid consisting of a mixture of multiple constituents by a 
reduced number of constituents (at least three), and a method allowing the opposite 
5 operation, a method referred to as delumping, i.e. a method allowing to obtain the 
detailed composition of the fluid from the description by a reduced number of 
constituents. 

The present invention thus refers to a modelling method : 

- allowing to estimate the properties or the behaviour of the liquid and/or vapour phases 
10 of mixtures of multiple components such as those encoimtered in situ in oil or gas 

reservoirs, or at the surface during production of these reservoirs, by means of a 
representation with a reduced number of pseudo-components, . 

- and allowing to predict, as a function of time, the detailed composition of fluids 
produced during production. 

15 The invention is particularly useful for accelerating calculations during simulation 

of the production of underground hydrocarbon reservoirs. 

Such a model allows reservoir engineers to reduce the time required for calculation 
of the behaviour simulations of reservoirs under production while keeping good 
modelling of the exchange mechanisms between hydrocarbon phases, and to establish 
20 detailed compositional profiles, which are necessary in particular for projected 
dimensioning and management of surface installations, such as separators, treating 
plants, transport lines, etc., therefore usefiil for surface and process engineering. 



BACKGROUND OF THE INVENTION 

Modelling of flows in an oil reservoir or in an underground storage is essentially 
based on the application to the previously gridded reservoir (or to a portion thereof) of 
the well-known Darcy's law describing the flow of fluids in porous media, of material 
5 balance laws in each volume unit, of thermodynamic relations governing the evolution 
of the phase properties of the fluids such as viscosity, density, on initial conditions, 
structural closure boimdary conditions and well conditions. 

The model known as "Black Oil", referred to hereafter as B.O., is one of the most 
commonly used models in petroleum simulation. It allows to describe a compressible 

10 three-dimensional and three-phase (water-oil-gas) flow. The petroleum effluents 
involved in this model are generally described by a water constituent, and two 
constituents for the reservoir fluid, the term constituent covering here the notion of 
component (as 1120 for water) and the notion of pseudo-component (grouping of 
components). The constituents involved in this model are three : a water constituent 

15 (W), a heavy hydrocarbon constituent (H) and a light hydrocarbon constituent (V). In a 
B.O. type model referred to as "strict", constituent (W) is present only in the water 
phase, constituent (H) is present only in the liquid hydrocarbon phase (referred to as oil 
or condensate), and constituent (V) is divided between the liquid and vapour 
hydrocarbon phases (gas phase). A B.O. model referred to as "extensive" differs from a 

20 "strict" B.O. model in that constituent (H) is divided between the liquid and vapour 
hydrocarbon phases. However, although the use of B.O. models is applicable to a large 
number of industrial cases, it is not advisable in a certain number of cases, in particular 
in the case of condensate gas reservoirs subjected to dry gas injection. 
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Another well-known simulation model, referred to as "compositional" model, is 
also used, wherein the hydrocarbon fluids are represented by a larger number of 
constituents, at least three, often more, the water constituent being present only in the 
aqueous phase. Modelling the flow of these more detailed fluids leads to calculating 
5 times that are all the longer as the number of constituents is great. 

In order to allow the modelling calculations to be carried out within a reasonable 
period of time, the fluids in place are described as consisting of a number of 
components or pseudo-components that is much more reduced than the real number of 
components. Switching from a detailed representation of the fluids to a representation 
10 with a smaller number of constituents is carried out by means of "lumping" or 
"pseudoisation" operations. In the description hereafter, unltess otherwise stated, the 
term "pseudoisation" is used for any method allowing to reduce the number of 
constituents. 

Various pseudoisation methods have already been proposed for selecting and 
15 defining the pseudo-components, and the engineer often has to find a compromise 
between precision and calculating time (and cost). For simulation of the produaion of 
condensate gas reservoirs subjected to dry gas injection, representations with about 6 or 
8 constituents are generally used, which leads to calculating times which are all the 
longer as it is often necessary to reduce the size of the grid cells to limit numerical 
20 errors and consequently to increase the number of grid cells. Considerable eflFort is 
devoted to the development of pseudoisation methods for use in the industry, which 
would allow to reduce even fiirther the number of constituents while modelling the 
behaviour of the fluids with precision, and making it possible to obtain detailed 
compositional information. The operations allowing to predict the reservoir simulation 
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results that would be obtained using finely detaiiled modelling (where the fluids are 
represented by a greater number of components) are known to the man skilled in the art 
as "delumping". 

Patent WO-00/37,898 describes a pseudoisation method applicable to compositional 
5 simulations, based on selection of a number of "dominant" base components equal to 
the number of pseudo-components desired at the end of the procedure. In this method, 
the non-dominant components are represented in all the pseudo-components, and a 
particular dominant component is represented in a single pseudo-component only. The 
mathematical transformation on which the lumping method is based allows, by inverse 
10 transformation, to obtain the detailed compositional information. The composition of a 
pseudo-component taken in particular can show negative molar fi-actions of base 
constituents, as in the example given in Table E7 of the patent mentioned by way of 
reference. It is understandable that such a representation loses a certain physical sense 
when a particular pseudo-component is considered individually. This may lead to 
15 robustness problems when, as it is the case in practice with gas injection, a "local" 
simulation result, for example in certain grid cells, shows the disappearance of one or 
more of the constituents used in the dynamic simulation. Besides, implementation of the 
invention is described as requiring many iterative calculations and a large storage space. 

Earlier publications describe pseudoisation methods also applicable to 
20 compositional simulations, wherein each pseudo-component is formed by grouping 
together several base constituents, a particular base constituent being represented in a 
single pseudo-component only. Lumping can be performed according to a selection of a 
priori set criteria such as those given in the aforementioned patent, or by means of an 
optimization procedure, for example as proposed by K. Liu in the paper "Reduce the 
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Number of Components for Compositional Reservoir Simulation", SPE 66363, 
presented at the SPE Reservoir Simulation Symposium, Houstoti, Texas, 11-14 
February 2001. 

A paper well-known to the man skilled in the art, vmtten by D.E. Kenyon and G. 

5 Alda Behie, "Third SPE Comparative Solution Project : Gas Cycling of Retrograde 
Condensate Reservoirs", SPE 12278, Journal of Petroleum Technology, August 1987, 
illustrates a situation that is not exceptional : a great disparity can be observed, for the 
same case study, in the compositional simulation results when the simulations are 
carried out with different fluid representations, and moreover with different simulation 

10 softwares. In this paper, the number of constituents ranges from 5 to 16, and the 
disparity of the results is in part due to the various selections of fluid compositional 
representations : the general tendency observed is that the smaller the number of 
constituents, the more the hydrocarbon liquid (referred to as oil or condensate) 
saturation in a particular grid cell can be underestimated, and the more the oil recovery 

15 at the surface is then overestimated. 

In order to collect the detailed compositional information during a compositional 
reservoir simulation, delumping methods such as those described in patent WO- 
99/42,937 and in the paper by C. Leibovici and J. Barker "A Method for Delumping the 
Results of a Compositional Reservoir Simulation" SPE 49068, presented at the SPE 
20 Annual Technical Conference and Exhibition New Orleans, 27-30 September 1998, can 
be used. The method allows to foresee the evolution of the detailed composition in time 
from calculations, in particular equilibrium coefficient calculations, carried out in a 
compositional type simulation of fluids described by a certain reduced number of 
pseudo-components, the number of components being at least three. 
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A Black Oil type representation can be considered to result from a particular 
pseudoisation operation providing two pseudo-components. The detailed composition of 
each pseudo-component, which is not useful for construction of the representation, is 
not known a priori, which is not a crippling obstacle in collecting detailed 
5 compositional information, by a delumping operation. Thus, patent FR-00/09,008 
describes a method allowing to foresee the evolution of the detailed composition in time 
from calculations carried out in a Black Oil type djoiamic simulation. 

The principle of the delumping stage of patents WO-99/42,937 and FR-00/09,008 is 
to determine, -from calculations carried out during simulation, with the lumped 

10 thermodynamic representation (compositional or BO), in each grid cell and at each time 
interval-, coeflBcient ADo and the n coefficients ADp (i.e. n+l coefficients, n being the 
number of parameters of the state equation) of a known general equation previously 
published in a paper by C.F. Leibovici, E.H. Stenby, K. Knudsen, "A Consistent 
Procedure for Pseudo-Component Delumping", Fluid Phase Equilibria, 1996, 117, 225- 

15 232 : 

n 

(1) Ln(Ki) = ADo + Z ADpI^ 

where IQ is the equilibrium constant of constituent i and the Tlpi are fixed parameters for 
characterizing constituent i in the state equation for a given thermodynamic 
representation. 

20 Once coefficient ADo and the n coefficients ADp determined, they are used to 

calculate the equilibrium constants of the constituents of the detailed thermodynamic 
representation (N^ components) by applying Equation (1) to the components with 



their own fixed characterization parameters in the detailed thermodynamic ^~ 
representation. 

One of the significant points of this method is that, in the delumping stage, it is not 
necessary to solve the Na equilibrium equations associated with the state equation 
5 (equations which express the equality of the fiigacities of each constituent in each 
phase) in the various time intervals of the dynamic flow simulation, which saves 
calculating time. 

The paper by W.H. Goldthorpe "Simulation of Gas Injection Processes in Gas- 
Condensate Reservoirs Using a Binary Pseudo-Component Representation", SPE 

10 19470, presented at the SPE Asia-Pacific Conference, Sydney, Australia, 13-15 
September 1989, illustrates the simulation results that can be obtained with an advanced 
Black Oil modelling in the case of production of a condensate gas reservoir by means of 
a gas injection process. The case taken as an example comes firom the aforementioned 
publication by D.E. Kenyon and G. Alda Behie. Considering the disparity of the results 

15 in this publication, the results obtained by W.H. Goldthorpe with a Black Oil 
representation, by comparison with the results of a simulation performed with a detailed 
representation, appear to be much more coherent, but it can be seen in Figure 8 of the 
paper that the solution is not satisfactory in the revaporization stage because the oil 
saturation in a particular grid cell (the same as in the reference paper) is very different 

20 from the saturation of the detailed compositional prediction, and seems to be truncated 
of negative values during eight simulated production years. 

The state of the prior art is thus such that there is no simple and robust lumping 
method available : 
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-. allowing to reduce to three the number of pseudo-components in compositional 
simulations, so as to obtain notably reduced calculating times, in particular for 
simulation of gas injection cases with revaporization effects, difficult to treat with a 
Black Oil representation, 

- guaranteeing a physical sense to the compositional simulation results and, 
consequently, to the associated delumping operation results. 

SUMMARY OF THE INVENTION 

The pseudoisation method according to the invention allows to estimate the 
properties or the behaviour of liquid and/or vapour hydrocarbon phases from data 
relative to a reference set consisting of hydrocarbon mixtures in a series of 
thermodynamic states resulting from the production conditions encountered, or 
expected as such, for underground hydrocarbon reservoirs. It comprises the following 
stages : 

- grouping each one of said hydrocarbon mixtures into at least three constituents (V, I, 
H), none of these constituents corresponding to a particular selection of base 
components or pseudo-components that would be used for a detailed compositional 
description of the fluids, considering that the gas phases resulting from the separation 
under surface conditions of each one of the hydrocarbon mixtures are mixtures from 
which third constituent (H) is excluded, and that the oil phases resulting from the 
separation under surface conditions of each one of the hydrocarbon mixtures are 
mixtures from which first constituent (V) is excluded, 

- determining by material balance the compositions of the separation products 
comprising, for the gaseous products, at least the first and the second constituent (V, I) 
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in variable proportions and, for the liquid products, at leaSt the second and the third 
constituent (I, H) in variable proportions, and 

- determining the at least three-constituent composition of each hydrocarbon mixture of 
the reference set by combination of the products of the separation thereof in proportion 
5 to the amounts of each separation product. 

According to an implementation mode, each one of the hydrocarbon mixtures is 
grouped into only three constituents (V, I, H), the gas phases resulting from said 
separation are mixtures in variable proportions of first constituent (V) and of second 
constituent (I), the oil phases resulting from said separation are mixtures in variable 
10 proportions of second constituent (I) and of third constituent (H), and the three- 
constituent composition is determined. 

The surface conditions are those encountered or expected during production of the 
reservoir, but they can be redefined depending on the context of the study carried out, 
and they can therefore be different from the surface conditions encountered or expected 
15 during production of the reservoir. 

According to an implementation detail, the material balance is a mass balance and a 
molar mass is assigned to each one of the three constituents (V, I, H) after quantitative 
analysis of the molar masses of the separation products of the reference set. 

According to an implementation variant, the data necessary for equihbrium 
20 calculation and for modelling the phase properties in the lumped representation are 
defined using the compositions of the phases in the lumped representation and known or 
estimated a priori data relative to at least the density and the viscosity of the oil and gas 
phases at equilibrium belonging to the reference set. • 
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When said data includes detailed compositional data of the phases previously 
represented by a "detailed" state equation, the parameters of a first (and possibly single) 
state equation of the lumped representation, used for modelling the phase properties, are 
defined using this compositional data. 

According to another preferred implementation variant, the parameters of a second 
state equation of the lumped representation, useful for equilibrium calculations, are 
adjusted in order to reproduce the equilibrium coefficients of the lumped representation. 

To perform this adjustment, the parameters are used for example per constituent of 
the lumped representation in a state equation useful for calculation of the phase 
properties. 

According to another implementation variant useful in particular for later 
delumping calculations, the equilibrium coefficients of the fluids are determined in a 
detailed compositional representation, from variables and/or parameters involved in the 
calculation of the phase properties, from the moment that the parameters useful for 
calculation of the phase properties in the lumped representation have been estimated so 
as to reproduce the parameters of the phases in the state equation of the detailed 
compositional description. 

The method can also comprise delumping stages for predicting as a fiinction of 
time, and in at least one thermodynamic zone, a detailed composition of a fluid 
contained in a hydrocarbon reservoir or produced by at least one well, these stages 
being for example as follows : 
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- representing the reservoir in form of a network of grid cells (m) wherein each one 
forms an elementary volume filled with fluids in form of one or more phases, with at 
least one non-aqueous phase, 

- defining, for each thermodynamic zone or range, the fluids by a detailed base 
5 representation, so as to determine the amount of each base constituent (i) in each 

hydrocarbon phase in each grid cell (m) at the time defined as initial for the delumping 
calculation, 

- per thermodynamic zone for which a lumped representation of the fluids is selected, 
determining a state equation constructed prior to dynamic reservoir simulation with the 

10 lumped representation, to reproduce the phase parameters, in the state equation of the 
detailed representation, of the hydrocarbon fluids along thermodynamic paths 
considered to be representative of those that will be followed by the hydrocarbon fluids 
during the gridded simulation, 

- carrying out, at a time interval t, a compositional simulation with a limited number of 
15 constituents ^ wherein the phase properties are calculated by a state equation, said 

simulation allowing to calculate at least in each grid cell (m) and at consecutive time 
intervals a pressure for a hydrocari>on phase, the temperature when it varies, the flow 

rates of the phases between grid cells, between perforated grid cells and well, and the 

J 

values of parameters and/or phase properties involved in the formal e?q)ression of the 
20 equilibrium coefficients of the detailed representation, and storing these various 
quantities, 

- estimating at the next time interval (t+1) the molar fi-action of each constituent i in the 
global detailed composition of the hydrocarbon fluid in grid cell (m) by material 
balance on grid cell (m). 
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- determining, using the quantities stored, at the same time interval (t+1) and in each : 
grid cell (m), the equilibrium coefficients of each constituent (i) in the detailed 
representation, 

- deterniining, in the same time interval (t+1), the vaporized fraction in each grid cell 
5 (m), and 

- estimating the detailed composition of each hydrocarbon phase, at the same time 
interval (t+1) and in each grid cell (m), in particular in form of the numbers of moles of 
each constituent i in the detailed representation in each phase at the same time interval 

(t-Hl). 

10 In other words, the essential points of the pseudoisation method can also be defined 

as follows : 

- Having data, at least density and viscosity data, available concerning oil and gas 
phases at equilibrium in situ for a certain number of thermodjaiamic states *e' under 
pressure and temperature conditions P*" and T*, 

15 - knowing directly or by means of previous calculations the molar masses of the oil 
phases resulting from the separation of the oil in situ, from the separation of the gas in 

situ, respectively denoted by (JV^o^)and {^m^^^^ the molar masses of the gas phases 

at the surface resulting from the separation of the oil in situ, from the separation of the 

gas in situ, respectively denoted by iMM^^^wA (iW^^), and the vapour molar 

20 . fractions (or gas phase molar fraction), G^^ for the separation of the oil in situ, 6^^ for 
the separation of the gas in situ, 

- determining the compositions of the phases "at the surface", and consequently the 
- equilibrium coefficients, in the lumped representation using constituents (V, I, H) with 
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the assumption that the gas phases at the surface can contam constituents (V) and (I), to 
the exclusion of constituent (H), and that the oil phases at the surface can contain 
constituents (I) and (H), to the exclusion of constituent (V), 

- determining the compositions of the phases "in situ" in the lumped representation 
5 using constituents (V, I, H) by recombination of the compositions of the phases at the 

surface using the vapour molar fractions e^^ , e^^ , these determinations providing the 
equilibrium coefHcients of constituents (V, I, H) under the "in situ" reservoir 
conditions, 

- using the a priori known data and the compositions of the phases in the lumped 
10 representation to define the input necessary for modelling the phase properties and for 

the equilibrium calculations in the lumped representation and, more particularly when 
the a priori data includes detailed compositional data of the phases previously 
represented by a state equation, using this additional information to define the 
parameters of a state equation of the lumped representation useful first for modelling the 
15 phase properties, and possibly for defining the parameters of a second state equation of 
the lumped representation useful for equilibrium calculations. 

Constituents V, I, H are respectively related to a volatile constituent, an 
intermediate constituent and a heavy constituent, without any of them corresponding to 
a particular selection of base components or pseudo-components that would serve for a 
20 detailed compositional description of the fluids. 

The essential idea implemented for delumping compositional simulations is to 
calculate the equilibrium coefficients of the fluids in the detailed compositional 
representation from variables and/or parameters involved in the phase properties 
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calculation from the moment that the parameters useful for calculation of the phase 
properties in the lumped representation have been estimated so as to reproduce the 
parameters of the phases in the state equation of the detailed compositional description. 

BRIEF DESCRIPTION OF THE FIGURES 

5 Other features and advantages of the method according to the invention will be 

clear from reading the description hereafter of a non limitative embodiment example, 
with reference to the accompanying drawings wherein : 

- Figure 1 diagrammatically shows the distinct separation of an oil phase and of a gas 
phase at equilibrium in the reservoir, 

10 - Figures 2 to 5 compare the results of a compositional simulation with a three pseudo- 
component representation (results referred to as ternary), obtained by applying the 
method according to the invention, with the results of a detailed compositional 
simulation with sixteen base constituents (results referred to as reference), 

- Figures 6 and 7 show, for the detailed and ternary representations, respectively the 
15 variation of the attraction term and of the covolume (expressed in dimensionless form) 

of the oil and gas phases as a fiinction of the pressure during the thermodjmamic 
simulation of a depressurization operation at constant volume at reservoir temperature, 

- Figures 8.1 to 8.16 on the one hand and Figures 9.1 to 9.16 on the other hand show, 
per constituent, respeaively the detailed composition of the gas and of the oil, during 

20 the same depressurization operation, as obtained from the reference representation, and 
from delumping of the three pseudo-component representation, 

- Figure 10 recapitulates the results illustrated in Figures 8.1 to 8.16 and 9.1 to 9.16 by 
showing the maximum absolute error obtained for the proportion of each constituent 
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during the depressurization from 237.4 to 69 bars as a fUnction of the proportion of 
constituent in the phase at the dew-point pressure of the initial fluid, 

" Figures 11.1 to 11.16 on the one hand and Figures 12.1 to 12.16 on the other hand 
show, per constituent, respectively the detailed composition of the gas and of the oil, 
during the thermodynamic simulation of a condensate vaporization operation by 
injection of a dry gas at a constant pressure of 169 bars, as obtained from the reference 
representation, and from delumping of the three pseudo-component representation, and 

- Figure 13 recapitulates the results illustrated in Figures 1 1.1 to 11.16 and 12.1 to 12.16 
by showing the maximum absolute error obtained for the proportion of each constituent 
for the vaporization operation as a ftmction of the proportion of constituent in the phase 
at the beginning of the vaporization operation. 

DETAILED DESCRIPTION 

Precanble 

Modelling the behaviour of fluid mixtures requires a definition mode for the 
constituents used to describe the composition of the fluids and comprises : 

- thermodynamic equilibrium calculations to know if, under given pressure and 
temperature conditions, the state of equilibrium of a mixture, characterized by a global 
composition, is a single-phase state, or if the mixture divides into different phases, the 
equilibrium calculation then allowing to know the proportions of each phase, and the 
composition of each phase. An equilibrium calculation on a mixture, whose global 
composition is described by the molar fractions Zi of each constituent (i), leads to 
determination of the vapour fraction 9 and determination of the composition of each 
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phase. The vapour fraction 9 is generally determined by solving the Rachford-Rice 
equation known to. specialists, which is applied for example in the following document : 

- Rachford, H.H.Jr and Rice, J.D. : "Procedure for Use of Electronic Digital Computers 
in Calculating Flash Vaporization Hydrocarbon Equilibrium", J. Pet. Technol, 1952, 
14, 19, 

/^i + (Ki-i)e " 

the compositions of the "oil" and "gas" phases, respectively described by the molar 
fractions Xi and yi , being determined by : 

Zj KjZj 

(3) Xi i + (Ki-l)e * l + (Ki-l)e 

A hydrocarbon phase is referred to as undersaturated at a given pressure and 
temperature T when it is not at equilibrium with another hydrocarbon phase. According 
to a convention often admitted by the man skilled in the art, an undersaturated 
hydrocarbon phase is referred to as "oil" (or "condensate") when its critical temperature 
is below temperature T, and it is referred to as "gas" when its critical temperature is 
above temperature T. When its critical temperature is equal to temperature T, the phase 
can be referred to, in principle indiscriminately, as "oil" or "gas", and it is therefore 
considered to be related to one or to the other "oil" or "gas" qualities. 

- phase properties calculations such as density, viscosity, 

- if need be, delumping calculations allowing to reproduce detailed compositional 
information when the fluids have been modelled in a non detailed way. 

The various stages of the modelling method proposed are presented hereafter, the 
presentation being structured in five subchapters : 
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- constituents definition mode, 

- phase properties, 

- thermodynamic equilibria, 
-delumping, 

5 - comments concerning the generalization of the approach. 
Constituents definition mode 

a) We consider a set "E", referred to as reference set, of non aqueous fluids, 
referred to as hydrocarbon mixtures even though they may contain certain components 
other than hydrocarbons, such as nitrogen, carbon dioxide, sulfur dioxide, likely to 
10 correspond to : 

- the undersaturated hydrocarbon fluid phases found, assumed or predicted in one or 
more thermodynamic zones of the reservoir considered, 

- the hydrocarbon fluid phases at equilibrium found, assumed or predicted in one or 
more thermodynamic zones of the reservoir considered, 

15 - the hycfrocarbon fluid phases corresponding to the hydrocarbon fluids injected (or 
assumed to be) in the reservoir. In the most frequent practical cases, the injection fluids 
are gases which mainly contain light and intermediate hydrocarbons. 

b) When the necessary data is available, a detailed representation with base 
components and/or pseudo-components is defined to characterize the fluid mixtures. 
20 This situation being very common, the rest of the presentation is placed in such a 
situation because it allows to progressively introduce the usual conventions usefiil for 
the description of the method. 

Number being typically above 10, each constituent is described by a certain 
number of parameters, including, its molar mass MMi, subscript i being the number 
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assigned to the constituent (i=l,..., Na). The molar fractions Xi and yi respectively 
describe the composition of the liquid ("oil") and vapour ("gas'*) hydrocarbon phases in 
the base representation. 

An element 'e' of set "E" is then described by : a pressure P% a temperature P*, and: 
5 - (2 X Na) molar fractions « .x;i>fr)^Or. .y^w,)when two hydrocarbon phases at 
equilibrium coexist, 

- or (Na) molar fractions (x^ ^.x'^^) or iy^,y^^) when a single undersaturated 
hydrocarbon (oil or gas) phase is present. 

The equations given hereafter are written for an element *e' where two liquid and 
10 vapour phases coexist. For an element *e' where a single undersaturated phase is 
present, the set of equations relative to the single phase represented simply has to be 
selected. 

c) We consider that we have, for each phase of the reference set "E", fluid data 
likely to correspond to the compositions of the found, assumed or predicted 

15 hydrocarbon fluid phases corresponding to the fluid mixtures produced by a separation 
chain, and the separation chain can be reduced to a single expansion or comprise one or 
more intermediate separation stages. A separation chain is here, as it is usual in 
reservoir engineering, considered to be defined by the pressure and temperature values 
of the separation stages. It is convenient to consider the existing or planned separation 

20 chain through which the fluid mixture is or is assumed to be produced. In some practical 
cases, different separation chains can be considered, for example a different separation 
chain per thermodynamic zone, or a separation chain change in a given thermodynamic 
zone. We simplify hereafter first the presentation by considering a single separation 
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chain identical for all the elements of set "E", then we consider the case of several 
separation chains. 

The separation products of an oil phase of composition (xf.jc*^) are generally a 
liquid phase (oil phase) and a vapour phase (gas phase) whose molar masses are 
5 assumed to be known and respectively denoted by MMJ^ and MM^^ , in the molar 
proportions assumed to be known of 0^ of vapour phase and (l-^o) of liquid phase. 

can take all the real values between 0.0 and 1.0, boundaries included. Subscript o (in 
bold type) reminds that these are separation products of the oil phase. Subscripts O and 
G (not in bold type) give the oil or gas quality of the separation product. 

10 The separation products of a gas phase of composition {yi y^rb) ^® generally a 

liquid phase (oil phase) and a vapour phase (gas phase) whose molar masses are 
assumed to be known and respectively denoted by MMgo and MMgg, in the molar 
proportions assumed to be known of Go of vapour phase and (l -6g) of liquid phase. Go 
can take all the real values between 0.0 and 1.0, boundaries included. Subscript G (in 

15 bold type) reminds that these are separation products of the gas phase. 

Figure 1 diagrammatically shows the distinct separation of an oil phase and of a gas 
phase at equilibrium in the reservoir in the general situation where the (simple or 
complex) separation of an initially single-phase mixture, liquid (oil) or gas, produces 
two phases, one gaseous, the other liquid. 

20 The following relations express the mass conservation during separation operations 

of respectively the oil and/or the gas from element 'e' : 
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NIM^ = (1 - Go) MMto + Go MMo^ 
MNfc = (1 - Gg) MMg^ + Og MMgo 

where MMo and MM^ are the molar masses of the oil and gas phases of an element 'e'. 

d) The pseudoisation method provided consists in describing each hydrocarbon 
phase of set "E" in at least three constituents (V, I, H). 

5 The pseudo-components are defined by considering a hypothesis denoted hereafter 

by "M", according to which the gas phases resulting from the separation under surface 
conditions of each one of the hydrocarbon mixtures are mixtures from which third 
constituent (H) is excluded, and that the oil phases resulting from the separation under 
surface conditions of each one of the hydrocarbon mixtures are mixtures from which 

10 first constituent (V) is excluded. 

Constituents (V), (T), (H) are respectively related to a volatile, intermediate and 
heavy constituent, without any one of them corresponding to a particular selection of 
base components or pseudo-components that would serve for a detailed compositional 
description of the fluids. The associated modelling is referred to hereafter as BRO, 
15 which is short for BROWN OIL, which evokes the colour of the oil that changes from 
black to brown with an increasing proportion of intermediate constituents. 

First of all, to facilitate the presentation, we consider, according to a restrictive 
interpretation of the master hypothesis "M", that the mixtures produced by the 
separation chain are : 



(4) 
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- a mixture in variable proportions of the volatile (V) and intermediate (I) pseudo- 
constituents for a gaseous product (the molar mass of constituent V being less than the 
molar mass of constituent I), and 

- a mixture in variable proportions of the intermediate (I) and heavy (H) pseudo- 
5 constituents for a liquid product (the molar mass of constituent I being less than the 

molar mass of constituent H). 

The restriction on hypothesis "M" will be lifted in paragraph tt). 
The following relations therefore apply : 



(5) f MMoS==x^iMMi + XohMMh 



MMo^ = yov MMv + yoiMMj 



MMgo Xgi MMi + Xgh MMh 

6 e c 

MMgg = ycv MMv + yea MMi 



Xoi + XoH = 1 

yov + yoi = 1 

X^I + X^H = 1 

yGv+yGi= 1 



For the ternary representation to retain a physical sense, it is necessary that : 

10 - molar masses MMv, MMi, MMh are defined positive, 

- molar fractions x^i,xoh, yov,yoi, x^i,34h, yGv,yGi range between 0.0 and 1.0, 
boundaries included. 

It results therefrom that the selection of molar masses MMv, MMi, MMh has to 
meet the following inequalities : 

15 • MMv Min {MMgg, MMqg} | 

Max {MMgg, MMog) I ^ MMi < Min {MMgo, MMoo} I 

|(e€E) KeeE) 
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• e -• e 

MMh ^ Max {MMgo, MN^bo) 



(eeE) 



In the previous equations, the braces denote a set of values on which a minimum 
(Min) or maximum (Max) value is sought. 

The inequalities being non strict, the absence of one of constituents (V, I) in a given 
5 gaseous separation product is possible, but it is not a general rule and, similarly, the 
absence of one of constituents (I, H) in a given liquid separation product is possible but 
it is not a general rule. 

e) The values of the molar masses of constituents (V), (I), (H) have to be selected 
while meeting the above conditions and they must all be different from one another. In 
10 the example given by way of illustration, we have set down, by analogy vsdth a Black 
Oil approach, without it being a Umitative embodiment example : 



MMv = Min {MMgg, MMog) I ; MMh = Max {MMgo, MMoo} 

I (eeE) 



(eeE) . 



The illustration example, widely used in the rest of the description hereafter under 
the designation of 'SPE3' case, corresponds to the first test case of the aforementioned 
15 publication by D.E. Kenyon and G. Alda Behie relating to the production of a 
condensate gas reservoir by gas cycling. 

f) Once the molar masses of constituents (V), (I) and (H) defined, the composition 
of the fluids from the separation chain considered at c) is deduced, by the hypotheses of 
paragraph d) and Equations (5), by : 
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(6) 



e — MMqo 



Xgv^O 



e _ MMi — MMgq 
\. y^^" MMi-MMv 



; e • MMh MMqo 



xoi 




- MN^ 


e 




-MMv 


> yoi 




MMv 




MMh- 


e 

MMgo 










MMgo 


-MMv 


; ygi 


MMi- 


-MMv 



XoH = 



MMqo-MMi 
MMh-MMi 



yoH = o 

e ^ MMgq — MMi 
MMh-MMi 

yGH = o 



g) The composition of the phases of element 'e' in the three-constituent 
representation is obtained using the following relations (where a notation in bold type is 
used to denote the molar fractions of the fluids in situ), which translate the molar 
conservation during separation operations respectively of the oil and/or the gas of from 
element 'e' : 



(7) 



Xk = (1 - Oo) XoK + Go yOK 
yK = ( 1 - Gg) + yOK 



K= V,I,H 



By introducing relations (6) in the set of equations (7), we obtain in the reservoir for 
an oil phase of element 'e' the following composition : 

W f e ^ e MMi-MMrc 

^^"^^MMi-MMv 

e e ' 

e_ e MMh - MMqo e MMoG - MMy 

xi-(l-eo) MMh-MMi MNt-MMv 



e 

^ Xh-(I-Bo) MMh-MMi 



10 and for a gas phase : 



24 



(9) 



e 

_ o MMi — MMgg 



MN4-MMv 



e e 
e_ e MMh - MMgo MMcq - MMy 

yi-(l-©G) Mi^-MMj MN^-MMv 



e 

_ _ e MMgo - MMi 



h) After stages a) to g), we thus have, for later use, minimum data for defining 
constituents (V), (I), (H), i.e. their molar masses. 



It can be noted that, for all of the elements 'e' corresponding to states of equilibrium 
between hydrocarbon phases, we have, from the sets of equations (8) and (9), the 
equilibrium coeflBcients of each constituent (V), (T), (H) : 



(10) 



e yv 9g MMi — MMgg 

Kv ~ e ~ e c 

xv MMj — MMoQ 

e e 

e MMh-MMco ^« MMgg ~ MMy 
e (^-^Qg) MMh-MMi '^^^ MMi-MMy 

^tCj — e — e c 

Xi e MMh - MMqo e MMqg - MMy 

^^"^^^ MMh-MMi ^ MMi-MMv 



e yn I-Og MMgo-MMi 
Xh I-OoMMto-MMi 



It can also be noted that the set of equations (6) defines the equilibrium coeflBcients 
for the outlet conditions, in pressure and temperature, of the separation chain considered 
10 in paragraph c) : 



(11) 



< 



ypy 

Koy= e =00 
Xov 



K^y = ¥^ = oo 
Xgv 



e yoi MMoa - MMy MMh - MMi 
e yct MM<^-MMvMNfa-MN<t 



XoH 



; k;h=^ = o 

Xgh 
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i) Set of equations (7), written for constituents (V) and (H) respectively, gives the 
following relations : 



(12) 



Xv — 6o yov 



y V = Og Ygv 



from which it follows 

(13) f e XV 

yov =-7 

Wo 



e Vv 

yGv=t^ 



then : 

(14) I 



e 



Koi — 



1- 



1-Go 



1- 



1- 



Set of equations (14) shows that, from the moment that the molar masses of 
constituents (V), (I), (H) are fixed, and that the compositions of the phases of an 
10 element *e' in the ternary representation have been established, the equilibrium 
coefficients of constituent (I) can be obtained from the sole knowledge of the molar 
fractions of constituents (V) and (H) in the oil and gas phases of element 'e' and of the 
vapour fractions 9o and 9g at the outlet of the separation chain for each phase. A first 



26 



implication is that, in case of multiple separation chains, it is possible to generate the 
equilibrium coefficient data of constituent (I) relative to a separation chain (referred to 
as CSJB hereafter) other than the chain used to obtain the oil and gas compositions of 
element 'e' (referred to as CS_A hereafter). This is diagrammatically shown in the 
5 following table : 



Detailed 
mixture 


Operation 


Data used 


Temaiy 
mixture 


Equiiibiiiun coefficients 
CS_B 

B 




CS_A 

^> 


e e e . 

{MNdoo, MMoG, 9o}cs^ 

. ^> 


e c c 

{XVj Xi, Xh}cs_a 


J Xvcs A 
« QoCS_3 

Xhcs a 




CS^B 
> 


e 

{0o}cS_B 





0 

1 — 0O CSJ 






^> 






(yi " Yna)^ 


CS_A 

— > 


e e e 

{MMgo, MMgg, 9g}cs^ 

^ > 


{yv, yi, 3^}cs^ 


J _yvcs^ 

^ yncs^ 




CS_B 

^> 


{Qg}cs_b 


_ ^> 


1 ~ ©GCS^ 






> 







The conditions required for such a procedure are the observance of the following 
inequalities : 

10 • yoics_B > 0 et xoics_B ^ 0 hence xvcs^a ^ 6ocs_b ^ 1 - xhcs_a 

yci cs_3 ^ 0 et xgi cs_b ^ 0 hence y v cs_a ^ Qg cs_b ^ 1 - yn cs^a 

It can be noted that the equilibrium coefficients of constituent (I) thus established 
are not formally equivalent to the coefficients that would be obtained by applying sets 
15 of equations (6) to (11) to separation chain CSJB, knowing molar masses {MM^, MMoo, 
MMoo, MMgo}cs_b, and these molar masses must satisfy, in relation to the molar masses 
of constituents (V), (I), (H), the inequalities described in paragraph d). 
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If nefed be, it is possible to review the selection of the molar masses of constituents 
(V), (I), (H) by taking into account the molar masses of the phases resulting from 
separation CS_B in the writing of the inequalities of paragraph d). 

j) By way of example, the non limitative procedure followed to obtain the BRO 
representation for the *SPE3' case presented as an illustration consisted in various 
stages : 

- generating with a thermodynamic simulator such as, for example, the simulator known 
as TVT PACKAGE', a set of states 'e/ by means of a depressurization operation at 
constant volume at the reservoir temperature (93.3*'C), from the fluid saturation pressure 
(237.4 bars) down to 69 bars, using a Peng-Robinson state equation, very commonly 
used within the petroleum context, for the detailed 16-constituent fluid whose initial 
composition and parameters are given in Table 1 hereafter : 
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Table 1 : SPE3 case - Reference description of the reservoir fluid 



Name 


Molar 


MM 


Tc 

op 


Pc 
DSi 


3 
ft^/lbm 


0) 






N2 


CI 


C02 


C2 


C3 






28.016 


-232.33 


492.84 


1.4433 


0.0350 




N2 


0.000 


# 




# 


# 






16.043 


-115.76 


672.98 


1.5858 


0.0130 




CI 


0.120 


0.000 


# 


# 








44.010 


88.07 


1069.7 


1.5057 


0.2250 




C02 


0.000 


0.15000 


0.000 


# 






.0869 


30.070 


90.34 


709.67 


2.3707 


0.1050 




C2 


0.120 


0.00000 


0.150 


0.000 


# 






44.097 


206.28 


61 7.28 


3.2037 


0.1520 




C3 


0.120 


0.00000 


0.150 


0.000 


0.000 






58.124 


274.96 


528.95 


4.2129 


0.1918 




IC4 


0.120 


0.00000 


0.150 


0.000 


0.000 






58.124 


305.64 


550.56 


4.1007 


0.2010 




C4 


0.120 


0.00000 


0.150 


0.000 


0.000 






71 151 


370.1 1 


483.41 


4.9337 


0.2060 




ICS 


0.120 


0.00000 


0.150 


0.000 


0.000 






72.151 


385.92 


489.36 


4.9817 


0.2520 




C5 


0.120 


0.00000 


0.150 


0.000 


0.000 




M R1 
.U1 01 




452.25 


468.47 


5.9268 


0.2809 




C6 


0.120 


0.03000 


0.150 


0.010 


0.010 




M AA 


100.205 


51 2.64 


396.68 


6.9200 


0.3520 




C7 


0.120 


0.03000 


0.150 


0.010 


0.010 


C8 


.0150 


114.232 


564.64 


362.02 


7.881 1 


0.3992 




08 


0.120 


0.03200 


0.150 


0,010 


0.010 


C9 


.0105 


128.259 


613.31 


331.99 


8.7621 


0.4439 




C9 


0.120 


0.03400 


0.150 


0.010 


0.010 


CIO 


.0073 


142.286 


655.07 


304.15 


9.6591 


0.4869 




C10 


0.120 


0.03600 


0.150 


0.010 


0.010 


C11 


.0049 


147.000 


699.91 


349.54 


10.172 


0.4770 




C11 


0.120 


0.04300 


0.150 


0.010 


0.010 


C12P 


.0138 


210.000 


797.52 


324.79 


14.000 


0.5266 




C12P 


0.120 


0.04408 


0.150 


0.010 


0.010 



In this table, the binary interaction parameters form a zero-trace symmetrical 
matrix; the boxes comprising a symbol # are filled by symmetry. The missing columns 
have to be filled on the one hand by the synmnietry property, on the other hand by zero 
binary interaction parameters. 

- At various depressurization stages distributed between the saturation pressure of the 
initial fluid and 69 bars, a constant-pressure (the pressure of the stage considered) 
vaporization operation of the oil phase of the state *ed' considered has been carried out 
by means of a gas whose composition is considered to be representative of the injection 
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gas during simulation (in the SPE3 case, the composition of the gas injected varies 
during simulation since it is the gas produced that is reinjected), each vaporization 
operation comprising several injected gas amount stages xmtil complete vaporization of 
the Uquid phase. A state 'Cdv' corresponds to each injected gas amount stage where the 
5 phases are saturated. 

- Still with the same simulator, the separation operations for the oil and gas phases of 
each state 'e/ and 'eav' forming set "E" have been simulated, the separation chain CS_A 
used comprising a first stage at 26.TC and 56.2 bars, a second stage of expansion of the 
liquid phase fi-om the first separation stage to 26.7°C and 4,5 bars, and a last expansion 
10 of the Uquid phase fi-om the second separation stage to 15.6''C and 1.01325 bars, the gas 
phases fi-om the various separations being mixed. Since the SPE3 case involves a 
modification, during simulation, of the pressure conditions of the first separation stage 
(the pressure changing from 56.2 bars to 21.7 bars), the separation operations of the oil 
and gas phases of each 'e' state have also been simulated with the other separation chain 
15 CS^B. 

It has to be noted that the sequence of the various simulation operations described 
above can be readily carried out with most of the thermodynamic simulation softwares 
available for industry. 

The required conditions described in paragraph d) found from the molar masses of 
20 the separation products of the oil and gas phases of each state of set "E" by chains 
CS_A and CS_B are in gram/mole : 

MMv< 21.592 ; 27.865 < MMi ^ 72.334 ; MMh^ 162.06 
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Thfe molar masses assigned to the three constituents in the BRO representation are 
given in Table 2. 

The three-constituent compositions of the phases of set "E" by the calculations 
described in paragraphs f) and g) are obtained by considering the separation products of 
5 chain CS_A. The composition of the initial fluid in the BRO representation is given in 
Table 2 hereunder : 



Table 2: SPE3 case - BRO description of the reservoir fluid 


Ckmstituent name 


MM 


Molar fraction 


. V 


21.592 


.856714 


I 


62.334 


.092358 


H 


162.06 


.050928 



Phase properties 

The Gibbs rule, which gives the number of degrees of freedom of a thermodynamic 
10 system from the number of constituents and the number of phases, applied to a three- 
constituent mixture, gives three degrees of freedom when two hydrocarbon phases are at 
equilibrium : the phase properties depend on the pressure, the temperature and a 
compositional variable. 

For each fluid phase for which the phase property values useful for modelling the 
15 flows (compressibility factor, viscosity,..) are known for given pressure and temperature 
conditions, these phase properties can be entered in form of tables which are a fimction 
of the pressure, the temperature, and an index of the composition in the three- 
constituent representation. This data can be used afterwards as input data in a simulator. 
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the simulator generally having internal interpolation or extrapolation methods for 
estimating the phase properties at intermediate points. 

k) Alternatively to the use of tables, it is possible to seek, for example by 
regression, the correlation or state equation parameters allowing to calculate, by 
thermodynamic zone, the properties of the phases as a function of the composition and 
the pressure and temperature conditions, these parameters being then introduced in the 
simulation model instead of the tables described above. 

A reliable state equation for the reservoir fluid described in detail in the base 
representation is frequently available. 

Thus, v/hen a state equation (and an associated correlation for viscosity calculation) 
defined for the detailed Nrt,-constituent representation allows to model the 
thermodynamic behaviour of the fluids during the pressure and temperature stages of a 
thermodynamic path followed in one or more thermodynamic zones, this state equation 
can be used to generate the reference" set of the compositions of fluids and the 
aforementioned properties of the corresponding phases, by adding thereto the 
parameters of the phases in the state equation, such as the covolume, the attraction term, 
the critical pressure,..., which can be considered to be specific properties of the phases. 
For each fluid composition obtained, the same state equation, or another one, can be 
used to generate the compositions and properties of the fluid mixtures produced by a 
separation chain, and thus allow application of the pseudoisation method described in 
paragraphs a) to i). 

When the parameters of the pseudo-components in a state equation, and in a 
correlation for viscosity calculation, are adjusted, for example by regression, to 
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reproduce the specific properties of the phases such as the covolume, the.attraction 
term, the critical temperature,..., obtained with the detailed representation, we obtain a 
state equation and a correlation that can be used for calculation of the phase properties 
such as the compressibility factor, the viscosity, but this first state equation, referred to 
hereafter in short form as EOS_PRO, is not necessarily valid for the equilibrium 
calculations notably if the observance of a condition of equality of the fugacities at 
equilibrium has not been taken into account. 

A non limitative procedure allowing to obtain the parameters of the state equation 
EOSJPRO of the BRO representation is the procedure followed for the 'SPE3' case 
presented by way of illustration. The state equation used, the Peng-Robinson equation, 
is reminded hereaft«- for better understanding of the approach, but the procedure can be 
readily applied to other state equations. 

1) The Peng-Robinson equation is written as follows for phase P : 

RT a(T)p 

(15) P=7j — r-- ^ J 

Vp-bp (Vp + bp) -2bp 

where R is the perfect gas constant, T the temperature, P the pressure, V the molar 
volume of phase P, b its covolume and a(T) its attraction term, the latter quantities being 
expressed by : 

(16) bp = 2:cibi 



(17) a(T)p = i:cj 



avec aij(T) = (1 - 5^ 



where Ci denotes the molar fi-action of constituent i in phase P, a molar fraction whose 
physical dependences are omitted to lighten the presentation, bi and a,<T) respectively 
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the covolume and the attraction term of constituent i, 5ij the parameter of binary 
attraction between constituents i and j, the 6ij forming a zero-trace symmetrical matrix. 

In dimensionless form, the Peng-Robinson equation is expressed as follows : 

3 2 ' 2 2 3 

(18) Zp --(l-Bp)Zp +(Ap-3Bp -2Bp)Zp-(ApBp-Bp -Bp) = 0 

5 where Z is the compressibility factor of phase P, B and A are the dimensionless forms 
of the covolume and the attraction term of phase P (the dependences of the variables 
with the pressure and temperature being omitted to lighten writing of the equations) : 

(19) Bp = bp^ 

P 

(20) Ap = ap 2 i 

R T 

10 m) The compressibility factor of phase P is only a function of the two dimensionless 

phase parameters A and B. It follows therefrom that estmiation of the parameters per 
constituent in the BRO representation (i.e. bv, bi, bn, av(T), ai(T), aaKT), 5vi, 5vh, 6ih) 
allowing to reproduce the values of a(T) and b, obtained for each phase at the various 
stages of the PVT simulations with the base representation (detailed with 16 

15 constituents), allows to reproduce the reference values of the compressibility factors of 
the phases. The regression procedure used to obtain the various parameters of the BRO 
representation is the minimization of objective functions of the form : 

2 g 2 

(21) 0(p) = Z ti4 [<|>g(P) - <4] + S ^ [q>o(p) - o^] 

e e 

where the observable quantity o per phase and per state e, obtained here with the 
20 detailed 16-constituent representation, is modelled by a functional cp parameterized by 
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the components p, of vector p, and where quantities ra are weights assigned to the 
various observables. 

The minimum of 0(p)is obtained by solving the system of equations, of same 
dimension as vector p, formed by the nullity of the partial derivatives of function 0(p) 
5 with respect to the various parameters p, : 

(22) ap^O(p) = 0 
i.e. : 

(23) £ tife [<Pg(P) - <4] q4(p) + Z tD^ [q)^(p) - o^] 9p,9^(P) = 0 

n) When seeking the values of parameters W (K=V,I,H) of the BRO representation, 
10 the objective function used in the SPE3 case is thus : 

2 g 2 

(24) 0(bv,bi,bH)= Z tDoj^EyK bK- ^ J tu^j^Z^K bic- boj 

The system to be solved is simply a linear system of three equations in three 
imknowns : 



(25) 



mw bv+ nivi bt+ nivH bn^Tv 
mp/ bv+ mn bi+ mm bn^Ti 
bv+ niffl bi+ Hte bn^TH 



15 the terms and rK of the second members of the equations being given by 



o) In order to seek parameters av(T), ai(T). anCT), 8vi, 6vh, 6m simply by solution of 
linear systems of three equations in three unknowns, equation (17) can be rewritten as 
the sum of two terms, the second term comprising the binary interaaions : 
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(26) a(T)p = ai(T)p^ + a2(T)p' 
with : 

(27) aiCr)p = Sci's/^ and a2(T)p = - L Z Q Ci aijOuaij(T) = 

i j > . 

In term a2(T)p, parameters aij(T) form a set of three independent parameters avi(T), 
avH(T), affl(T). 

Using terms ai(T), obtained during the PVT simulations with the detailed 16- 
constituent representation, as observable quantities o per phase and per state, the 
unknowns being quantities ^ay{T),^aj(T),^aji(T) , the problem to be solved is actually 
similar to the problem described to obtain the covolume parameters of the constituents 
of the BRO representation. Similarly, parameters aij(T) are obtained by solving the 
system of three equations in three unknowns formed from observable quantities a2(T). 

p) The binary interaction parameters can then be obtained by : 

^ avTm avH(T) ^ olth(T) 

^^^^ ^-4^ 

q) Some simulators accept, as parameters of the constituents for the state equation, 
direct introduction of the covolumes bi of the constituents and of the matrix of 
parameters a4j(T) involved in the calculation of the attraction term (see equation (17)), 
and this matrix has to be defined if need be for different temperatures. 

If the case treated reqiures a temperature dependence (variable temperature in the 
reservoir),the previous procedure can be adapted, for example by generating set "E" by 
PVT operations at various reservoir temperatures, considering states e at various 
temperatures to obtain the covolumes per constituent (which have to be the same for the 
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different temperatures), considering the states e obtained at a given temperature to 
generate the corresponding matrix of parameters aij(TO and by repeating the procedure 
for as many temperatures as necessary. 

r) When the simulator does not accept direct introduction of the covolume per 
5 constituent bi and of the matrices of parameters ay(T*) : 

- the covolume per constituent is generally obtained from the introduction of parameters 
TCi, Pci, optionally Qbi (by default all equal to o!077796) by : 

(29) bi^RQbi^ 

- quantities ai(T) are generally obtained by : 

10 (30) V^C^5=^/aOi|^l + Q[l-^J^] with Q^^Cn^!" 



Quantity aOi is generally obtained from the introduction of parameters TCi, Pci, 
optionally Qai (by default all equal to 0.457235) by : 



(31) a0i = R2Qai||- 

In polynomial Q of degree dp, parameters cOi are the acentric factors of the 
15 constituents generally considered to be constant, and the degree of the polynomial and 
constants Cn can vary from one simulator to the next. 

s) When the critical pressures, temperatures and volumes per constituent are 
required, the critical pressures, temperatures and volumes of the phases of set "E", 
defined by mixture laws identical to those used for the covolumes, can be used as 
20 observables : 



2 
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(32) Pcp^LciPCi ; - Tcp = ZciTCi ; Vcp = i:ciVCi 

i i ' 

Parameters Pcv, Pci, Pch, Tcv, Tpi, Tch and Vcv, Vci, Vch respectively can thus be 
obtained by successive solution of three systems of equations in three unknovms totally 
similar to those generated for calculation of the covolumes per constituent, developed 
above. 

t) Parameters Clb^ (K=V,I,H) are then obtained from equation (29) and parameters 
bK, Pck and Tck as determined in paragraphs n) and s), by : 

PCKbK 

(33) ObK=R^ 

u) In equation (30), in the most frequent case where the reservoir temperature is 
considered to be constant, the temperature dependence of quantities 
4ai^.,laJ(T),^{a^W) can be disregarded by selecting a zero value for the acentric 
factors of constituents V, I, H, hence : 

(34) Vi;;(T5=>/aO^,VK,K=V,I,H 



and parameters (K=V,I,H) are then obtained from parameters (T) , Pck and Tck 
as determined in paragraphs o) and s), by : 

(35) ^^K= r^tck" 

Another possibility, used in the SPE3 case, consists in assigning to parameters Qbk 
(K=V,I,H) a value by defauh. In this case, parameters V^o]7 are determined with the 
PCk and Tck as determined in paragraph s), and parameters ^ are then obtained by 
rearranging equation (30) : 
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(36) Ck = 




T_ 



1 

then acentric factors ©k are sought as solutions to the three polynomial equations 
(K=V.I,H). 

(37) icn(i>K -Qi = 0 

v) When, for the case considered, the reservoir temperature is variable, equation 
(30) can be rewritten as follows : 

(38) VaKCr) = pK + YKnK(D 



vvrith the quantities : 'fK(T') = l-^-^ calculated using the critical temperatures of the 
constituents as determined in paragraph s). 

We can then seek quantities pK and Yk ifix = 4^ Tk = 4^ ) as independent 



parameters by linear regression of the pairs of values iiKiT'),4°KiT"X^^^ latter values 
being determined in paragraphs o) and q). 

Parameters QaK and are then obtained by : 

<39, Oa..^ and t = g 

The acentric factors sought are solutions to the three polynomial equations 
(K=V,I,H) : 

(40) ZCna)/-^=0 



n=0 
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w) Volume correction factors are often used to adjust the phase densities. When 
volume correction factors per constituent are required, it is possible to use as observable 
quantities the volume corrections of the phases of set "E", defined by mixture laws 
identical to those used for the covolumes : 

5 (41) CVp = ZciCvi 

i . • 

Parameters Cvy, Cvi, Cvh can then be obtained by solving a system of three 
equations in three unknowns totally similar to the system generated for calculation of 
the covolumes per constituent developed in paragraph n). 

x) In the SPE3 case, the hydrocarbon phase viscosity calculations are carried out 
10 with a correlation widely used in the field of petroleum engineering, the Lohrenz, John, 
Bray, B.G. and Clark, C.R. correlation : "Calculating Viscosities of Reservoir Fluids 
fi-om their Compositions", Journal of Petroleum Technology, 1964, pp. 1171-1176. The 
phase parameters explicitly used in this correlation are the critical pressures, 
temperatures, volumes of the phases. Determination of the critical parameters (pressure, 
15 temperature, volume) of constituents V, I, H as described in paragraph s) has allowed 
suitable modelling of the phase viscosities with the BRO representation. The procedure 
can be readily applied to other correlations. 

y) The guiding line followed throughout paragraphs n) to w) consists in seeking the 
state equation parameters EOS_PRO usefiil for calculation, under the reservoir 
20 conditions, of the phase densities (via the compressibility factors) and of the viscosities 
by simple solutions (non iterative) of systems of linear equations. 
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Its efficiency in reproducing the covolumes and attraction terms of the phases is 
illustrated, for the fluids of the SPE3 case, by Figures 6 and 7 for a depressurization 
operation at constant volume at the reservoir temperature (93.33°C) from the dew-point 
pressure bars of the initial reservoir fluid (237.4 bars) to 69 bars. Each figure shows the 
5 quantities calculated from the state equation of the detailed reference representation and 
the quantities calculated from the state equation EOS_PRO of the BRO model, these 
quantities practically merging with one another. 

This same line has been successfiiUy used for treating other cases that are not 
described here. It has been observed that such a line can be widely adapted to variable 
10 temperature contexts. Besides, the necessary processings pose no particular coding 
problems and have been easily automated. 



The parameters of the EOS_PRO used for the SPE3 case illustrated are given in 
Table 3 hereunder : 



Table 3: SPE3 case - Parameters of the EOS_PRO of the BRO representation 




Name 


Tc 

op 


Pc 

psi 


Vc 

tf/lbm 




Qb 


Cv 
tf/lbm 






v 


I 


H 


V 


-61 .212 


666.02 


1.869 


-0.02043 


0.079345 


-0.0686 




V 


0.0 


0.058447 


0.002065 


1 


368.292 


581 .21 


4.389 


0.31255 


0.082877 


4).0621 




1 


# 


0.000000 


0.001604 


H 


663.138 


263.46 


10.99 


0.21233 


0.065003 


0.2664 




H 


# 


# 


0.000000 



15 z) In the SPE3 case, we have chosen to use for the EOS__PRO a state equation of 

the same type as the equation used for the base representation. This is not obligatory. It 
is possible to use a state equation of another type, from the moment that this other state 
equation uses phase parameters identifiable with the phase parameters used by the 
reference state equation. 
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'More generally, 'the parameters required for the BRO representation can b.e obtained 
by various techniques known as regression, inversion, optimization,..., that can be found 
in an extensive literary output. They can for example be sought by rheans of one or 
more iterative solutions of systems of equations formed by partial derivatives, in 
5 relation to subsets or to the set of parameters of the BRO representation, of an objective 
function of the type of equation (21) or of another type. 

We could even try to use a formalism of this type to optimize one or more of the 
molar masses of constituents V, I, H by using it (or them) explicitly - by replacement in 
the objective function of the phase compositions by their expressions given by 
10 equations (8) and (9) -, and by selecting observable quantities sensitive to their value(s). 
Such an approach could be considered in a situation, which has not been encoimtered so 
far, where the guiding line followed throughout paragraphs n) to w) leads to a 
detennination of the BRO representation parameters that does not allow accurate 
reproduction of the observable quantities. 

15 Finally, for reasons known per se, one may be led to give a priori values to some of 

the parameters of the EOS_PRO equation and then to carry out estimation calculations 
only for a limited number of parameters. 

aa)The first expression generally expected from the simulation results are volume 
results. The observables which then have to be reproduced are the molar volumes of the 
20 phases after separation. For the SPE3 case, we used two distinct correlations, one for 
calculation of the molar volumes of the gas phases at the separation outlet, the other for 
molar calculation of the oil phases at the separation outlet, the separation outlet 
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conditions being the standard conditions used in the. profession, i.e. an atmosphere at 
15.56X (60°F), 

Many correlations or state equations are available for estimating the molar volume 
of a gas phase, notably under conditions close to atmospheric conditions. For 

5 information, it is the P.M. Dranchuk and J.H. Abou-Kassem correlation (presented in 
the paper : "Calculating Z Factors for Natural Gases using Equations of State", JCPT, 
July-September 1975) that has been used for the SPE3 case where the parameters usefiil 
for calculation of the molar volume of the gas are the critical pressures and temperatures 
of constituents (V) and (I), parameters which have been subjected to an iterative 

10 regression so as to reproduce the molar volumes of the gas phases obtained by the PVT 
separation simulations with the 16-constituent representation. 

The variation, as a function of the molar fraction of constituent I, of the molar 
volume of the oil phases under conditions close to atmospheric conditions has been 
found, in the SPE3 case, practically linear, more generally, for cases that are not 
15 described here, representable by a polynomial of low degree (degree two to three). For 
the SPE3 case, it is a linear dependence that has been selected, whose parameters are 
obtained by simple linear regression. 

The parameters required for calculation of the molar volumes under separation 
conditions were obtained for the two separation chains CS_A and CS_B with the same 
20 method. 
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Thermodynamic equilibria 

bb)The equilibrium calculations can be carried out either from equilibrium 
coefficient data, or from the expression of the equality of the fugacities of all the 
constituents in each phase, these two options being respectively presented hereafter. 

cc) According to Gibbs' rule, for calculation of the equilibria in the reservoir, the 
equilibrium coefficients of pseudo-components (V), (I), (H) as calculated by the set of 
equations (10) can be introduced in tables as functions of the pressure, temperature, and 
of an index of the composition in the pseudo-representation. This data can then be used 
as input data in a simulator, the simulator generally having internal interpolation or 
extrapolation methods for estimating the equilibrium coefficients at intermediate points. 
It is also possible, prior to introducing the data in a simulator, to add points to the tables 
using interpolation or extrapolation methods known in the art, or to use, instead of the 
tables, correlations constructed to reproduce the equilibrium coefficients. 

A dependence of the equilibrium coefficients with the pressure, the temperature and 
the composition index is not always necessary. Thus, in many application cases, one 
may be led to leave out the dependence with respect to a compositional index. These 
cases, where the number of degrees of freedom of the equilibrium coefficients would be 
reduced, remain within the scope of the method provided. 

dd) For calculation of the equilibria at the end of a separation chain, the equilibrium 
coefficients of pseudo-components (V) and (H) are respectively infinite and zero, 
whatever the pressure and temperature conditions at the separation chain outlet and 
whatever the composition of the mixture subjected to separation. If need be, in a 
simulator that may not accept the concept of infinite or zero equilibrium coefficients 
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(this concept being often implemented by forcing the absence of the constituent 
respectively in the oil phase and in the gas phase), it is possible to use a very .great finite 
value and a very low finite value, both positive, instead of infinite or zero equilibrium 
coefficients. For the data relative to the equilibrium coefficients of constituent (I) 

5 necessary for equilibrium calculation, the set of equations (11) gives a different 
definition of the equilibrium coefficients of constituent (I) according to whether the 
separation of a 'gas' phase fi-om the reservoir or the separation of an 'oil' phase from 
the reservoir is considered. The input data, in principle necessary for equilibrium 
calculations, is thus, per phase, a table of equilibrium coefficients of constituent (I) as a 

10 fiinction of a compositional index of the phase subjected to separation, of the pressure 
and of the temperature at the outlet of the separation chain considered. The variation 
ranges of the compositions of the oil and gas phases being generally distinct, it may be 
convenient to "merge" into a single table the data useful for the oil and gas phases, 
notably if the simulator does not accept introduction of distinct data for each phase. 

15 This data can be used as input data in a simulator, the simulator generally having 

internal interpolation or extrapolation methods for estimating the equilibrium 
coeflRcients at intermediate points. It is also possible, prior to introducing the data in a 
simulator, to add points to the tables using interpolation or extrapolation methods 
known in the art, or to use, instead of the tables, correlations constructed to reproduce 

20 the equilibrium coefficients. 

ee) If the EOS_PRO state equation is constructed to reproduce the phase properties, 
it is not necessarily valid for equilibrium calculations, i.e. to allow reproduction of the 
compositions of the phases of set "E" obtained in paragraph g). As above, the Peng- 
Robinson equation is used hereafter as presentation support. 
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For this 'State equation^ the equilibrium coefficients, according to the condition «f * 
equality of the fugacities of each constituent in each phase at equilibrium, are obtained 
by: 

N 

(42) Lii(Ki) = Do + D,bi+;^|||7:^^(^yj-|^J9)as(T) 

j-i 

5 Avith the dimensionless phase parameters defined according to equations (19) and (20) 
and : 

(43) f , ^Zp.(l.V2)Bp- | 

^P + (l->/2)Bp>/ 

Pr Zo-l _ Zo-l AqLq _ AgLq" 

^ ^^"RtL Bo Bo 2^2 Bo^ iViBqI 

and the compressibility factors of the phases calculated from equation (18). 

For a case where the reservoir temperature is constant, the problem consists in 
determining the nine independent parameters forming a vector p : K, bi, bn, aw, an, ann, 
10 avi, avH, am, for example by minimizing an objective function constructed using as 
observables the equilibrium coefficients of set "E" as calculated by equations (10). 
The objective fimction according to equation (21) can then be adapted as : 

^ 2 

(44) 0(p) = Z E [<pkp) - oL] 

where functions tpf^ (p) are the exponentials of the second members of equations (42). 
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Alternatively, in equation (44), we could use as observables Oj^ the logarithms of 
the equilibrium coefficients of set "E" and as functions <p^ (p) the second members of 
equation (42). 

By following K.E. Starling's procedure in the paper : "A New Approach for 
Determining Equation-of-State Parameters Using Phase Equilibria Data", SPE 1481, 
SPEJ, December 1966, we can select weights crj^ as the inverses of the values of 
observables O^r . 

The nullity of the partial derivatives of function O(^) with respect to the nine 
parameters of interest provides a system of nine non-linear equations in nine unknowns, 
that can be solved iteratively only, many methods suited to solution of non-linear 
systems being discussed in the literature. 

An iterative calculation of this type requires a first approximation of the parameters 
vector It has been observed, notably for the SPE3 case, that it is appropriate to profit 
from the parameters obtained for the EOS_PRO to start the iterative calculations. 

More generally, notably to reproduce equilibrium coefficients at various 
temperatures, we can consider another set or subset of parameters for vector p and 
other parameter seek techniques. 

ff) By way of illustration, for the SPE3 case, vector p consists of the six parameters 
Pck and Qbk, the values of the other parameters being taken firom equation EOS_PRO, 
The procedure used for the SPE3 case is a procedure known in the art that is being 
improved. The set of parameters EOS^EQ used for the SPE3 case is given in Table 4 
hereafter: - 
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Table4:SPE3cas 


»e - Parameters of theEOS.EQof the BRO representation 




Name 


Tc 

op 


Pc 

psi 








5kk' 


V 


I 


H 


V 


-61.212 


848.47 


-0.02043 


0.089345 




V 


0.0 


0.058447 


0.002065 




368.292 


274.23 


0.31255 


0.070000 




1 


# 


0.000000 


0.001604 


H 


663.138 


309.90 

1 


0.21233 


0.065003 




H 


# 


# 


0.000000 



gg) In the SPE3 case, the data introduced for the surface equilibrium calculations 
for each separation chain CS_A and CS_B are the equilibrium coefiBcients of 
constituent (I) in form of a single table for the separation of the oil and gas phases, a 
5 table obtained as described in paragraph dd) with, as the single abscissa, the index of the 
composition of the phase before separation. 

hh) Figures 2 to 5 compare the results, over the fifteen simulated years of 
production of the SPE3 case, of a compositional simulation with three pseudo- 
components BRO (results referred to as ^ternary'), where the compositional 
10 representation was obtained by applying the method according to the invention, with 
those of a detailed compositional simulation with sixteen base constituents (results 
referred to as 'reference'). 

The reference compositional simulation uses the same Peng-Robinson state 
equation for the equilibrium calculations and for the phase property calculations, and for 
15 all the thermodynamic zones. The three-constituent BRO simulation uses, for 
calculation in the reservoir grids, state equation (EOS_EQ) for the equilibrium 
calculations and state equation (EOSJPRO) for the phase property calculations. For the 
surface condition calculations, the BRO simulation uses equilibrium coefficients 
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introduced in form of tables ' and sim^ile correlations for calculation of the phase 
volumes. 

The results shown in Figure 2, which represent the evolution of the pressure at the 
bottom of a producing well, and in Figure 5, which represent the oil (or condensate) 
5 saturation evolution in the deepest perforated grid in the producing well, are among the 
most revealing results of the right calculation of the equilibria and phase properties in 
the in situ conditions. 

It can be observed that the BRO model provides a pressure solution that practically 
merges with the reference solution (with 16 constituents) and an oil saturation evolution 

10 solution in good agreement with the reference solution over the total simulated time. 
The grid used in Figure 5 is the same one as the grid used in the reference paper by D.E. 
Kenyon and G. Alda Behie and in the paper by W.H. Goldthorpe mentioned in the prior 
art. It can be observed that the results of Figure 5, judged by comparison with the results 
obtained by W.H. Goldlhorpe with a BO modelling, show a much better quality of 

15 reproduction of the oil saturation evolution in the reservoir, whereas the procedure used 
to obtain the parameters of EOS_EQ is improvable. 

ii) The results shown in Figures 3 and 4 for the oil flow rate and cumulative oil 
production throughout the simulated time show that the BRO modelling provides a 
solution that is very close to the reference solution. 

20 The calculating time reduction provided by the BRO modelling is considerable, the 

acceleration factor being 46. 
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Delumping - 

jj) We first rqjresent the reservoir in form of a network of grid cells (m), each one 
forming an elementary volume filled with fluid(s) in form of one or more phases, with 
at least one non-aqueous phase. The non-aqueous phases are still referred to as 
5 hydrocarbon phases even though they may contain certain components other than 
hydrocarbons, such as nitrogen, carbon dioxide, sulfur dioxide. 

kk)We define, for each thermodynamic zone or range, the fluid by a detailed 
representation with components and/or pseudo-components. It can be noted that it is 
possible to treat the cases which require several thermodynamic representations, for 
10 example if several local thermodynamic paths can be distinguished during modelling 
(zones produced only by depressurization and zones subjected to gas injection can be 
distinguished). Several variation zones or ranges for the thermodynamic or 
compositional quantities, often referred to by specialists as thermodynamic zones, can 
thus be defined and used. 

15 11) Per thermodynamic zone for which we select a lumped representation of the 

fluids, we determine a state equation EOS JPRO constructed, prior to dynamic reservoir 
simulation with the lumped representation, to reproduce the phase parameters, in the 
state equation of the detailed representation, of the hydrocarbon fluids along 
thermodynamic paths considered to be representative of the paths that will be followed 

20 by the hydrocarbon fluids during the gridded simulation. 

mm) We carry out, in a manner known in the art, a compositional simulation 
with a limited number of constituents where the phase properties are calculated by a 
state equation EOS_PRO, said simulation allowing to calculate at least in each grid cell 
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(m), and* at consecutive time intervals (t, t+1, etc.), a pressure in a hydrocarbon phase 
{pi), the temperature (r^)(if it varies), the molar flow rates per injection or 
production phase (5^) and (SlJ, for each pair of grid ceUs (m,h), the molar flow rates 
of the liquid (m„^0 and vapour (u^) phases, the vapour fraction (^^). These quantities 
5 are stored, as well as the parameters in state equation EOS_PRO for the liquid and 
vapour hydrocarbon phases (typically the covolumes and the attraction terms) involved 
in the e5q)ression of the equilibrium coefficients. In order to avoid recalculating them, 
the compressibility factors of the phases (Z^) and (Z^) obtained by solution of state 
equation EOS_PRO can also be stored. 

10 For a Peng-Robinson state equation, while keeping the notations introduced by 

equations (16), (17), (26) and (27), the phase parameters are denoted by (bol) and (bom) 
for the covolumes of the phases, (acta) and (aol), - or (ancL), (aiom), (a20m), i&iom) for 
the attraction terms of the phases. 

The various parameters useful for equilibrium coefficients calculation, when 
15 mentioned in a general way, are denoted hereafter by (E^l ) , subscript P being for each 
oil and gas phase, and subscript tt a reference number for a particular parameter. 

nn) We estimate at the time interval t+1 the molar fraction of each constituent i in 
the global detailed composition (z,^') of the hydrocarbon fluid in grid cell (m), knowing 
the number of moles of each constitu«it i in the detailed representation in each phase at 
20 the time interval t, respectively iNol,)foT the oU phase and (Ngl^) for the gas phase 
from the equations hereafta* : 
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'' ■ .'Nib Nib t ' t t"' ' ' 

(45) Nol = SNOil ; Ngl=ZNsl ; Nm = No„. + Na„ 

i i 

t NOil t Ngil t NOil + Nal 

^^^^ ^"^Noi ' ^^^Ng: ' ^' n: 

ZilNl- At(yil' + SoL)- At Z (Jil' "gmh ^il' 

h6T(m) 

t+1 

(47) Zi^ -til : [ 

(48) =nI- At(S^ + ScL)- At E (Ugrfi + Uolh) 

heJCm) 

5 Equations (45) and (46) respectively allow to calculate, at time interval (t) and in 

grid cell (m), the number of moles per hydrocarbon phase and the global number 
Nl&om knowledge of the number of moles per constituent and per phase, the molar 
fractions of the constituents per phase and the global molar fractions. 

Equation (48) translates the total hydrocarbon molar balance on grid cell (m), by 
10 taking account of the material exchanges, during time interval At, with all the grid cells 
(h) adjoining (m) which form set J(m). 

In equation (47), the writing of terms {y^*)cmd (x^ - wherein m'=m for a flow rate 
from grid cell (m) to grid cell (h) or in the well, and m'=h for a flow rate from grid cell 
(h) to grid cell (m), and m' corresponding to the fluid injected in the case of injection 
15 wells, S being then negative- implicitly implies the use of a simple upstream scheme for 
the compositional flows. A more general writing of these terms is iy\^)cmd (x,^), 
where and describe the compositions of the Uquid and gas phases, obtained in 
a manner known in the art, flowing between grid cells (m) and (h). 
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10 



oo) We deteraiine, at time interval (t+1) and in each grid cell (m), the equilibrium 
coefficients {Kl^^)of constituent i in the detailed representation using a formulation 
explicitiy involving parameters Em», by introducing therein the variable parameters 
(jE^^J^^) resulting from the simulation with the lumped representation and, if need be, 
using some of the compositional variables described in the previous paragraph, 
calculated at time interval (t). Thus, for a Peng-Robinson state equation, the equilibrium 
coefficients given by equation (42) can be rewritten in the following form : 



t+i t+i _ t+i _ t-hi ^ / ^t+i P 

(49) Ln(Ki„,) = Dom bi + Dsm Vai(T„) +;^^2^t 



2 _t+l 2 J-im 
m 



with : 
(50) ^ 



✓ t+l trt\ 
Dom "^^nlr^S t+l I 

\Zom — Bomy 



t+l 
Lpm 



^ P=0,G 



t+l t+l 

Apm Lq 



.+ (1-V2)Bp1v 

t+l r M-i „ t+l 

^ t+l Pm_ Zom - 1 Zom - 1 

'^'"=^L Bo^' " < "2>/2Bor^"2>/2B, 

t+l r t+l 
t+l -Pm I Lpm 

D2m 2 tfTa ^ t+l 



t+l t+l 



t+l 2 — Z n _ t+l 2 I 



t+l ^ tfl 

raiGm ~ t+l aiom 
Bom 



Equations (50) involve the dimensiohless phase parameters B and A previously 
defined by equations (19) and (20). 

In equation (49), term //^^ represents the contribution of the binary interactions 
between constituents and this term is zero when the binary interaction coefficients of the 
15 detailed representation are all zero. In such a case, the equilibrium coefficients can be 
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calculated directly from parameters {Ej^^ ) from the state equation EOS JPRO of the 
dynamic simulation using the lumped representation. 

In the opposite case, a rigorous expression for term II^ is : 

t+i t+i 

Bom j=i i^Gm j-i 

5 This term thus involves the molar fractions, in the detailed representation, of the 

constituents in the hydrocarbon phases at tune interval (t+1) , not estimated 

yet. In order to avoid a heavy implicit resolution, it is possible, in equation (51), to 
substitute the molar fractions x/^\>'/^^ for the molar fractions at time interval (t) from 
equations (46). 

10 It is with such a non limitative approximation that the delumping results : 

- illustrated in Figs. 8.1 to 8.16, Figs. 9.1 to 9.16 and recapitulated in Fig.lO for a 
depressurization operation of the fluid of the SPE3 case, 

- illustrated in Figs. 11.1 to 11.16, Figs. 12.1 to 12.16 and recapitulated in Fig.13 for a 
vaporization operation at a pressure of 169 bars of a condensate (the condensate 

15 obtained at 169 bars during the previous depressurization operation), 

were obtained. 

In Figs. 8.1 to 8,16 and 9.1 to 9.16 on the one hand. Figs. 11.1 to 11.16 and 12.1 to 
12.16 on the other hand, the proportion of each constituent in the phase considered is 
expressed in molar percentage of the phase and, in each figure, the scale is expanded to 
20 be suited to the variation of the proportion of constituent in the phase. 
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The abscissa in Figs. 11.1 to 11.16 and Figs. 12.1 to 12.16 is the molar fraction of 
vapour during the vaporization operation, which increases from 0 (in the initial state, 
only the condensate phase, or oil, is present) to 1 when all of the condensate is 
vaporized. 

Fig. 10 and Fig. 13 show that the maximum absolute error obtained for the 
proportion of each constituent (still expressed in molar percentage) remains below 
0. 1 % on average during the depressurization and vaporization operations. 

We thus observe that the compositions resulting from the delumping operation are 
consistent with the detailed reference compositions of the fluids, although the 
approximation used is not among the most sophisticated ones. 

Another possible, still non limitative approximation consists in estimating the molar 
fractions x.'^\y'l^ from the equilibrium coefficients Kl^ that can be calculated from 
the molar fractions given by equations (46) by : 
(52) r t^i ^ zl^ 

""'"^ i-h(K^-i)er 

t+l _ kil Zim (Kim - 1) 



In general, the binary interaction coefficients are small in relation to 1, and an 
iterative calculation, for example by successive substitutions, is not necessary but 
remains possible. 

pp) We determine, at each time interval (t+l), the vaporized fraction (^l^^)in each 
grid cell (m), either from the results of the dynamic simulation with the lumped 



55 



representation or, if a higher precision is desired, by solving the Rachford^Rice equation . 
from the molar fractions of each component i in the global detailed composition 
(r/^' ) of the hydrocarbon fluid in grid cell (m) at time interval (t+1). 



(53) Z *^'",,^l. =0 
w l+(Ki„-l)9„ 

The numbers of moles per phase in each grid cell at time interval (t+1) are given 

by: 

tH t+1 t+1 ^ t+1 tfl 

(54) No„ =N„. (l-e„) ; Ngi„,=N^0„, 

qq) We estimate the detailed composition of each hydrocarbon phase, at each time 
interval (t+1) and in each grid cell (m), using the following relations : 

t+1 t+1 

(55) x'^ = Ix fertile oil phase 

l+(Kim-l)Om 

t+1 t+1 ^t+i 

(56) yZ = ^^" "t+r for liie gas phase 

l + (Kin-l)e„ 

rr) Finally, we obtain the numbers of moles of each constituent i in the detailed 
representation in each phase at time interval (t+1), respectively (iVb,^' ) for the oil phase 
and {Ngl^^) by equations (46) modified into : 

t+1 t+1 t+1 w-l t+1 t+1 

(57) NOim = NOm Xim \ Ngim = Ngm Yun 

ss) Paragraphs nn) to rr) detail the content of the operations for shifting from time 
interval t to time interval t+1, assuming that the numbers of moles of each constituent i 
in the detailed representation in each phase at time interval t in each grid cell of the 
simulation with the lumped representation, respectively denoted by(M>;„)for the oil 
phase and and (iVg^) ^r the gas phase, are known. To start calculation, these 
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quantities just have to be known at a given time or time interval^ which can be in 
particular the initial time. This can be obtained by various means known in the art, one 
means consisting in an initialization calculation of a compositional simulation 
(calculation of the initial state only) using the detailed compositional representation, the 
5 inputs other than compositional being the same as for the compositional simulation with 
the lumped representation. 

We then know how to describe the evolution of the detailed composition in each 
grid cell during the production process modelled in a compositional model using a 
lumped representation by means of a memory of the parameters of the state equation of 
10 the detailed representation kept in the phase properties calculation method. 

Some calculation stages may appear unnecessary because redundant. Some of these 
redundancies have the advantage of allowing to reduce the storage space required for 
data storage, but they can obviously be avoided if the storage space considerations are 
not a priority. 

15 In the case where EOS_PRO is valid for calculation of the equilibrium coeflScients 

(Ki, m and t+l), delumping of the results of a reservoir simulation carried out with the 
method according to the invention can be performed by applying the method described 
in the aforementioned patent WO/42,93 7. 

Comments concerning the generalization of the approach 

20 tt) The details provided in the description of the method use a molar formalism; a 

mass formalism could also be used. 
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uu) It can be noted that the. line of reasoning based on the master hypothesis "M" 
described in paragraph d), applied to the set constituents of the base representation, a set 
referred to hereafter as "B", in the previous paragraphs, could have been alternatively 
applied to a subset of the constituents of the base representation, referred to as "R" 
hereafter. The complementary set of 'TR." in "B" being denoted by "S", equations (5) 
would then be rewritten (with P=0,G) as : 



(58) 



ieS 



MMpo = XpiMMi + 34hMM:h + Z) ?4MMi 

i 

e e e e 

MMTO = ypvMMv + ypiMMi + X ypiMMi 



e c ^ e 

3^ + XpH=l-2- ^ 



; ypv + ypi= 1-2- 



Such an approach can be interesting when, for the case treated, it is necessary to 
keep explicit certain base constituents in the lumped representation. Thus, in a practical 
case where the reservoir conditions would not allow to disregard the dissolution of some 
constituents (such as carbon dioxide) in water, these constituents could form set "S". 



The simple variable changes given hereafter allow to bring equations (5) back to the 
situation where the set of constituents of the base representation has been considered : 



(59) 



ieS 



MMre-S y^MMi 



MMpoL = 



MMpgL = 



l-Zx^i 



ieS 

i-Zy^i 



< 



XpT 



XpH 



1-2 x^ l-2xpi 



c ypv 
ypv|R = — 



i€S 

i-Zy« 



m. 



i-Zy« 



58 

Equations (6) being solved with the new variables, an inverse transformation of / 
variables then allows to carry out one after the other the calculations of paragraph g) 
and the following paragraphs. 

The lumped representation then comprises constituents (V), (I), (H) and the base 
5 constituents of set "S". In each phase of set "E", the molar fractions of the base 
constituents forming set "S" are the same in the lumped representation as in the base 
representation. ■ 

w) The compositions of the phases of set "E" in the lumped representation being 
thus obtained, the broad lines of the stages described above can be followed to model 

10 the phase properties and the thermodynamic equilibria. From the moment that a state 
equation is sought for modelling with the lumped representation the phase properties 
and/or the thermodynamic equilibria, the indication is to retain for the base constituents 
forming set "S" the parameters of these constituents in the state equation of the detailed 
description. The only parameters to be estimated are then the parameters of constituents 

15 (V), (I), (H) in the state equation, and if need be the parameters of the binary 
interactions of these constituents with the constituents of set "S". When the latter are 
considered to be zero, the guiding line followed throughout paragraphs n) to w) can be 
followed, which consists in seeking the state equation parameters EOSJPRO useful for 
calculation, under the reservoir conditionis, of the phase densities (via the 

20 compressibility factors) and of. the viscosities by simple solutions (non iterative) of 
systems of linear equations. 



